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We discuss how contributions to the order 0(ma : ) orthopositronium 
, decay rate can be separated into two categories, one due to relativis- 

ed ' tic momenta and calculable in terms of QED scattering amplitudes, the 

other due to low momenta and calculable in the simpler framework of 
■ a low-energy effective theory (NRQED). We report new results for all 

O |' low-momentum contributions, and give a formula relating the remaining 

I , contributions to conventional (on-shell) QED scattering amplitudes. 

Ph. 

Despite its many successes, quantum electrodynamics has yet to account 
'-^ ' fully for the decay rate of orthopositronium. The theoretical expression for this 

> '• rate is §, § 
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where the coefficient B has not yet been computed. The measured rate is 
Aexp = 7.0482(16)/is _1 fiL The difference between this value and the known 
part of A t h, 

\ exp - X th (B = 0) = 99(16) x 10-V" 1 (2) 

is surprisingly large: the coefficient B in X t h would have to be about 250 to 
bring theory and experiment into agreement. Such a large coefficient would be 

'Present address: Dipartimento di fisica teorica deirUniversita' di Torino, Via P. Giuria 1, 
10125 Torino, Italy. 
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unusual, but is by no means impossible, particularly given the big 0(a/n) cor- 
rection. A complete calculation of B is essential before any realistic assessment 
of the situation is possible. 

A calculation of these coefficients using traditional bound-state methods is 
very complicated. This is because each term in a traditional expansion has 
contributions from both nonrelativistic and relativistic momenta. This means 
that approximations only valid for small p or only for large p cannot be readily 
employed to simplify the analysis. In this paper, we outline a new and simpler 
procedure for computing B. Our analysis is based upon a rigorous nonrel- 
ativistic reformulation of QED called nonrelativistic quantum electrodynamics 
(NRQED) @. Using this effective field theory, we are able to separate the 0(a 2 ) 
corrections into three parts. Two of these involve soft, nonrelativistic momenta 
and probe the bound-state nature of the system. The other part involves hard, 
relativistic momenta and is therefore largely insensitive to the details of binding. 
We have calculated the nonrelativistic contribution and present the results here. 
We also show how to extract the relativistic contribution from a calculation of 
ordinary scattering amplitudes; no bound-state physics is required in this part 
of the calculation. 

Our positronium results also have implications for quarkonium decays. These 
will be discussed in another paper 

The lowest-order decay rate of orthopositronium (O-Ps, n = J = S = 1) is 
given by (we now use natural units, with c = Ti = 1) 

r (O-Ps^3 7 ) = |*(0)| 2 <7 (0) 

2tt 2 - 18 R 
= m P a 

= 7.2112 fjs- 1 , (3) 

Here ^(O) is the ground state Schrodinger-Coulomb wavefunction evaluated at 
r = 0, and <7o(p) is proportional to the lowest-order decay rate of a free electron 
and positron in an S-state: 

<7 (p) = -L Im Mt° } (ee -> 3 7 -> ee) , (4) 

where p is the magnitude of the electron center-of-mass momentum. 
To order 0(a 2 ), there are three sources of corrections: 

1) RADIATIVE CORRECTIONS TO <t (0): These corrections renormalize 
the lowest-order decay rate: 

STi = SZ r (O-Ps -» 37). (5) 

The "renormalization" constant SZ is defined to be the part of the radiative 
corrections coming from relativistic loop momenta. Thus it is insensitive to the 
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binding energy of the atom, and can be related to the infrared-finite part of 
the radiatively corrected annihilation rate a for a free electron and positron at 
p = 0: 

<x(0) - cr IR (0) = (l + SZ)a (0), (6) 

where <tir(0) is the infrared part of <r(0). Using a photon mass to regulate the 
infrared (threshold) singularities^ 

^# = ^ + (21n2 + l)« 2 !f 
cr o (0) A \ 2 

, . a 2 m„ 1 A , o, 

- 2 x 10.282(3) + -a 2 In — + 0(a 3 ). (7) 

7r A 3 m e 

These infrared subtractions remove all contributions due to nonrelativistic loop 
momenta. Thus 5Z has a A-independent expansion in powers of a/ir: 

5Z = C i- + c 2 f-) 2 + .... (8) 

The only 0(a) corrections to the decay rate come from 5Z, and thus, from Eq.|jJ 
c\ = —10.282(3). The second-order coefficient has not yet been calculated. The 
large size of ci suggests that c-i might also be large. Indeed, a small (gauge 
invariant) subset of the diagrams contributing to c-i has recently been evaluated 
and found to contribute 28.8(2) to C2 @. 

2) MOMENTUM DEPENDENCE OF a (p) : Near threshold, the decay 
rate for a free electron and positron in an S-wave has the form 

5-0 (p) = a (0) + Sa ( P ) + O(^a ), (9) 



where 5<7q (p) can be computed to be Q : 



p 2 I9n 2 - 132 . 
5a {p) = —— — — cr 0). 10 

?7lg 12(7T 2 ™ 9) 

This correction term shifts the decay rate in 0{a 2 ). Naively, one might expect 
a contribution 



d 3 p 
(2^)3 



2 *(°) / 77TZT3^o(p)*(p). (11) 



However, this integral includes contributions already present in the relativistic 
radiative corrections discussed above (Eq.||). This becomes apparent if the 
Schrodinger equation is used to rewrite the wavefunction: 



d 3 q 

E _ Ps -p 2 /m e J J2nf 



*(P) = ~ / 77^3 ^(P " q) HQ) (12) 



1 The precise definition of ctir(0) depends upon the regulator used in the NRQED part of 
the analysis; see refs [S. 
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where Vc is the Coulomb potential, Vc(p — q) = — e 2 /(p — q) 2 , and #0-Ps = 
— ma 2 / 4 is the binding energy. Then our contribution to the orthopositronium 
decay rate becomes 



2 * (,)1 I (S^ 8vo(p)m = 2*(0) 1 ''''' 



J 44 S Sa C (E . Fs ,p)^(p) (13) 



where 



Sao(E Q . Ps ,p) . /(ga^odP + ql) ^,^^ ^) (14) 

is a one-loop radiative correction to the basic decay. 

To avoid double-counting the relativistic radiative corrections, we subtract 
the threshold value of Sac, 

5a c (E,p)-*6ac(E,p)-5a c (0,0). (15) 

Thus the final contribution to the decay rate from Sa (p) is 

6T 2 = 2*(0) J -^(6a c (E _ Ps ,p)-6ac(0,0))y(p). (16) 

Since S&oip) ocp 2 , this result simplifies to 

E Q _ Ps 19^-132 

« 5 x 10-V 1 - (17) 

Note that our original expression for this contribution (Eq.Ju]]) is ultravio- 
let divergent. The divergence comes from the high-p region, which is precisely 
the region that must be removed so as to avoid double counting the relativis- 
tic radiative correction. Consequently, our subtraction regulates the integral. 
Such infinities are common when nonrelativistic expansions like that for ao(p) 
(Eq.|J) are employed. However, they can always be removed by appropriate 
counterterms. The structure of these counterterms is readily determined using 
NRQED. The use of this theory puts calculations such as ours on a rigorous 
foundation, permitting systematic improvement of the final results (a rigorous 
derivation of the NRQED expansion to second order in p/m and a, including 
all counterterms, will be derived in another paper ||). 

3) CORRECTIONS TO *(0) : The wavefunction at the origin in the lowest- 
order rate (Eq.j^j) is modified by the various relativistic corrections that must 
be added to the Coulomb- Schrodinger theory. Using first-order perturbation 
theory, the correction is 

sm = E *4o)™f> (is) 

^O-Ps E 0-P S ~ E n 
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where the sum is over all eigenstates of the unperturbed theory, and 

, 3e 2 e 2 e 2 (q • <r x ) (q • <r 2 ) 



(T 1 ■ CT 2 



Am 2 ^m 2 4m 2 q 
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-3^1n(<z 2 /m 2 ). (19) 

The potential SV contains all 0(v 2 /c 2 ) corrections to the Schrodinger theory 
as well as the leading 0(av 2 / c 2 ) correction, given by the ln(g 2 /m 2 ) term. This 
last term gives the leading contribution to the Lamb shift in positronium. 
The corrections to the wavefunction will shift the decay rate by 

(5r 3 (0-Ps -> 37) = 2|£#(0)$(0)| ct-q(p). (20) 

As in the case of 6a(p), the ultraviolet divergences encountered in evaluating 
(^(O) are systematically removed by using NRQED. The remaining contribu- 
tions all come from nonrelativistic loop momenta. We find 

59(0) = [ia 2 lna + 0.58a 2 - ^-a 3 (lna) 2 ] *(0) (21) 
which implies 

ST 3 = [ia 2 lna+ 1.16a 2 - J-a 3 (lna) 2 ] r (O-Ps -> 3 7 ). (22) 

The logarithmic terms agree with the literature || |j; the other correction is 
new. 

By combining the 0(a 2 To) contributions from all three of our corrections, 
we obtain the final result 

= 9.56 x 10^/is" 1 +c 2 0.39 x 10" 4 /xs _1 . (23) 

The coefficient c 2 is specified by the threshold rate for ee — > 37 (Eqs.rol 0, raj). 
The known part of c 2 M contributes 11 x 10~ 4 /is _1 to the rate. 

In this paper, we have outlined a new procedure for analyzing the O-Ps 
decay rate at O(o; 2 ro). We have computed all corrections in this order that 
depend in detail on bound-state physics. The only remaining contribution can 
be extracted from a calculation of the annihilation rate for a free electron and 
positron - no bound-state physics is required. Our corrections combined with 
the known part of the c 2 correction account for 20% of the difference between 
theory and experiment. The unexpectedly large size of these corrections makes 



5 



it plausible that the full calculation of C2 will bring theory into agreement with 
experiment. 
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